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Abstract 

There are several known data structures that answer distance queries between two arbi- 
trary vertices in a planar graph. The tradeoff is among preprocessing time, storage space 
and query time. In this paper we present three data structures that answer such queries, 
y—i ' each with its own advantage over previous data structures. The first one improves the query 

^-H , time of data structures of linear space. The second improves the preprocessing time of data 

' structures with a space bound of 0{n'^^'^) or higher while matching the best known query 

^ ^ , time. The third data structure improves the query time for a similar range of space bounds, 

' at the expense of a longer preprocessing time. The techniques that we use include mod- 

. ifying the parameters of planar graph decompositions, combining the different advantages 

' of existing data structures, and using the Monge property for finding minimum elements of 

matrices. 

(N . 

C/Q^. 1 Introduction 

^ • There are several known data structures that answer distance queries in planar graphs. We 

O , survey them below. All of these data structures use the following basic idea. They split 

I the graph into pieces, where each piece is connected to the rest of the graph only through 

CN ■ its boundary vertices. Then, every path can go from one piece to another only through these 

, boundary vertices. The different data structures find different efficient ways to store or compute 

^ I the distance between two boundary vertices or between a boundary vertex and a non-boundary 

OO ' vertex. 

. Frederickson jl6j gave the first data structures that can answer distance queries in planar 

I graphs fast. He gave a data structure of linear size with 0(n log n) preprocessing time that can 

f"*) ' find the shortest path tree rooted at any vertex in 0{n) time, where n is the number of vertices 

^— I ■ in the graph. This leads also to an 0{n?) solution to the all-pairs shortest-paths problem, and 

^ ', implies a distance query data structure of size 0{n?) with 0(1) query time. Feuerstein and 

Marchetti-Spaccamela [15] modified the data structure of [16] and showed how to decrease the 
5^ ' time of a distance query by increasing the preprocessing time. They do not provide an analysis 

of their data structure in terms of preprocessing time, storage space, and query time, but they 
show the total running time of k queries, which is O(nfe-l-nlogn), 0{n^^^k^^^), 0{n^^^), 0{ny/k) 
for k < -y/n, ^/n <k<n, n<k< n^^^, n^^^ < k < , respectively. This solution actually 
consists of three different data structures for the three cases k < -^/n, -y/n < k < n and n < k, 
where the data structure for the first case is the one of [16] . 

Henzinger, Klein, Rao and Subramanian [19] gave an 0(n) algorithm for the single-source 
shortest path problem. This implies a trivial distance query data structure, which uses the 
algorithm, and takes 0(n) space and query time. 

Djidjev [9] gave three data structures. We will use the specific section number in |9] - §3, 
§4, or §5, to refer to each one of them. The first one [U (§3)] works for S G [n^/^,n^] and has 
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size 0{S), 0{S) preprocessing time, and O^ii^/S) query time. The same data structure was 
also presented by Arikati et al. [3]. This data structure is similar to the two data structures of 
[TB] , but takes advantage of the algorithm of [TO] to get a better preprocessing time. The second 
[g (§4)] works for S € [n,n^/^] and has size 0{S), 0{nVS) preprocessing time, and 0{n'^/S) 
query time. The third data structure [9l (§5)] works for S G [ri^/'^,n^], has size 0{S), 0{n\fS) 
preprocessing time, and 0{n\og{n/^fS)/^fS) query time0 Chen and Xu [8] presented a data 
structure with the same time and space boundsH 

Fakcharoenphol and Rao HD gave a data structure with 0(n log n) space and 0{n\og'^ n) 
preprocessing time and 0{\/n\o^ n) query time. Klein [22] improved the preprocessing time 
of the data structure to 0(n log^ n). 

Cabello [5] presented a data structure that uses 0{S) space and can be constructed in 0{S) 
time for S € [n^/^ log^^^ n, n^], and answers a query in 0{n\o^^'^{n) / y/S) time. This data 
structure can answer k queries in a total of 0(n^/^ log^^^ n + /c^/^n^/^ log n) time. If the queries 
are known in advance, the algorithm of [S] avoids storing the entire structure, and uses only 
0{n + k) space. 

Other data structures exist for outerplanar graphs (some use the fact that these graphs 
have treewidth 2) and planar graphs with small face-on- vertex covering [U [6l [TOl W\\ Wl\ I18j . 
for planar graphs with bounded edge lengths [13^ I24j. and for approximate distance queries in 
planar graphs [S E [SH [22l [30] . 

In this paper we present three new data structures for the problem: 

• SectionO A data structure with 0{n) space, 0{n log n) preprocessing time, and ©(n^/^"*"^) 
query time, for a constant e > 0. This data structure has the best known query time 
achievable with linear space. This also improves the total running time for answering 
k distance queries when k is 0(n^/^~'^) and a;(logn), and if we limit ourselves to data 
structures with 0(n+fc) space then the upper bound on the range of k grows to 0(n^/^~^). 
The data structure is based on the data structure of [2], by replacing the recursive 
decomposition of [H] with recursive r -decomposition as in [16]. As the data structure 
of [H], our data structure also generalizes to graphs embedded in a surface of bounded 
genus. 

• Section HI For S e [n^/^,n^], a data structure with 0{S) space, 0(51ogn) preprocessing 
time, and 0{nlog{n/\fS)/\/^) query time. This data structure matches the storage space 
and query time of [8l[9l (§5)], which is the best query time for this range of storage space, 
with a better preprocessing time for 5 = o(n^/log^n). The data structure is obtained 
by combining a preprocessing algorithm similar to [5] with a data structure similar to |9l 
(§5)]. 

• Section[5l For S G [n^/^, n^], a data structure with 0{S) space, requiring 0{{S'^/'^ / y/n) log^ n) 
preprocessing time, and 0{n/y/S) query time. This data structure has the fastest query 
time for the same range of storage space as the previous data structure, but with a longer 
preprocessing time. The fast query time is obtained using an efficient minimum search in 

a Monge matrix. 

The different data structures are summarized in Table [TJ 

^Djidjev [5] (§5)] presents the result for S £ [n^^^ ,ri?^'^] with 0{n\og{n) / \fS) query time, however the same 
data structure works for a larger range of S, and the running time is actually 0{nlog{n/y/S)/y/S). 

^Chen and Xu [S] do not bound the time and space of the data structure in terms of S, the bounds here are 
derived by setting r = jS in the bounds that appear below Lemma 28 (page 477) of [Sj; the bounds stated in 
[8] depend on the minimum number p of faces required to cover all vertices of the graph (p is called face-on-vertex 
covering), these bounds are obtained using hammock decomposition [T7], which can be applied to any planar 
distance data structure. 
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Reference 


Storage space 0{S) 


Query time 


Preprocessing time 


mm (§3)] 


Elil (§5)] 


o — n 

S e [7l4/3 logl/3 71, 7l2] 
Se [71^/3, 7l2] 


0(7l2/5) 

Oin\og'/^{n)/^/S) 
0{n\ogin/VS)/VS) 


U(n ) 
0{S) 
0{S) 
0{nVS) 


Section E] 
Section [5] 


S G [71^/3^ j^2] 

r /I 9i 

5 e [71^/3, 7l2] 


0{n\og{n/^/S)/^/S) 
0{n/VS) 


0(5" log 7i) 
0{{S^i^/^) log'^n) 


m (§4)] 
im 


e [71, 7l3/2] 

S = n log 71 

S* = 71 
71 


Oiy^log" n) 
0{n) 
0{n) 


0{ny/S) 
0{n log^ 7i) 
0(71 log 71) 


Section [3] 


= 71 


0(7ll/2+^) 


0{n log 71) 



Table 1: Comparison of distance query data structures for planar graphs. Time bounds are 
expressed as a function of the storage space. The data structures are ordered by decreasing 
storage space and then by decreasing query time. 

2 Preliminaries 

We consider a directed simple planar graph G{V, E). We let n = \ V{G)\^ and by Euler's formula 
|£'(G)| = 0{n). We assume that G is given with a fixed planar embedding, in other words it is 
a plane graph. Without loss of generality we assume that G is a triangulated, bounded degree 
graph; this assumption is required by algorithms that we use and are described in Sect. I2.11 and 
12.21 We assume that G is connected, since we can handle each connected component separately. 

Every edge in E{G) has a non-negative length. The length of a path is the sum of lengths 
of all of its edges. The distance from a vertex u to a vertex v is the minimum length of a path 
from u to V. With additional 0(n log^ n/ log log ti) preprocessing time we can allow negative 
edge lengths as weh, see [HI [231 EH] for details. 

Let F,H be subgraphs of G. We write dH{u,v) to denote the distance from u to f in H. 
The graph FDH is the subgraph induced by E{F) n E{H). For short we denote \H\ = \ V{H)\. 

2.1 Decomposition 

A decomposition of a planar graph G is a set of subgraphs of G, such that each edge is in exactly 
one subgraph and each vertex of G is in at least one subgraph. Each of the subgraphs which 
define the decomposition is called a piece. 

A vertex u is a boundary vertex of a piece if u G V{B) and v is incident to some edge not 
in E[B). The set of all boundary vertices of B is the boundary of B, denoted by dB. A hole 
is a face of B (including the external face) that is not a face of G. For a hole H we denote by 
H also the subgraph of G inside H. A boundary walk of i? is a facial walk of B around a hole 
H. For a piece B with hole H we denote dB[H] = dB D V{H). A vertex of V{B) \ dB is an 
internal vertex of B. 

All distance query data structures mentioned in the introduction decompose the planar 
graph. They take advantage of the fact that a path can go from one piece to another only 
through boundary vertices. 

A recursive decomposition [13] is obtained by starting with G itself being the only piece in 
level of the decomposition. At each level, we split each piece B with \B\ vertices and \dB\ 
boundary vertices that has more than one edge into two pieces, each with at most 2|i?|/3 vertices 
and at most 2\dB\/3 + 0{y^\B\) boundary vertices. We require that the boundary vertices of 
a piece B are also boundary vertices of the subpieces of B. Each piece B in the decomposition 
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has 0{y^\B\) boundary vertices. 

An r -decomposition [16] is a decomposition of the graph into 0{n/r) pieces, each of size at 
most r with 0{^/r) boundary vertices. 

Fakcharoenphol and Rao [H] showed how to find a recursive decomposition of G, such that 
each piece is connected and has at most a constant number of holes. They use these two 
properties for their distance algorithm. The construction of the decomposition takes O(nlogn) 
time using 0(n log n) space, and is done by recursively applying the separator algorithm of 
Miller |26j . Frederickson showed how to find an r-decomposition in 0(n log n) time and 
0(n) space by recursively applying the separator algorithm of Lipton and Tarjan [25]. Thus, an 
r-decomposition is a limited type of recursive decomposition where we stop the recursion earlier 
(when we get to pieces of size r), and do not store all the levels of the recursion (we store only 
the leaves). Cabello [S] combined the two constructions of |14l I16| (using [2B] instead of |25| ) 
and constructed an r-decomposition with the properties that the number of holes per piece is 
bounded by a constant, and that each piece is connected. 

In Sect. [3] we use a combination of recursive decomposition and r-decomposition - we de- 
compose the graph recursively, but we decompose each piece into 0{n/r) pieces instead of 2. 
In Sect, m we use r-decomposition. In Sect. [5] we use r-decomposition as well, there we take 
advantage of the fact that the construction of an r-decomposition is the same as of a recursive 
decomposition, which was stopped earlier. 

2.2 The Dense Distance Graph 

Fakcharoenphol and Rao [TJ] define the dense distance graph of a recursive decomposition. For 
each piece B in the recursive decomposition they add a piece to the dense distance graph that 
contains the vertices of dB and for every u,v € dB an edge from uio v whose length is dsiu, v). 
The multiple-source shortest paths algorithm of Klein [22] finds k distances where the sources 
of all of them are on the same face in 0{{k + n) log re) time. Therefore, using [22 it takes 
0(|i3| log \ B\) time to find the part of the dense distance graph that corresponds to a piece B 
(recall that \dB\ = 0{y^\B\) and B has a constant number of holes). It thus takes O(relog^re) 
time to construct the dense distance graph over all pieces of the recursive decomposition. 

Every single edge defines a piece in the base of the recursion, so it is clear that the distance 
from u to w in the dense distance graph is the same as the distance between these two ver- 
tices in the original graph. Fakcharoenphol and Rao noticed that in order to find the distance 
from u to V we do not have to search the entire dense distance graph, but that it suffices to 
consider only edges that correspond to shortest paths between boundary vertices in a limited 
number of pieces. The pieces are these containing either u or v, and their siblings in the recur- 
sive decomposition. There are O(logre) such pieces with a total of 0{^/n) boundary vertices. 
Fakcharoenphol and Rao gave an implementation of Dijkstra's algorithm that runs over a sub- 
graph of the dense distance graph with q vertices, defined by a partial set of the pieces in the 
recursive decomposition, in 0{qlog'^ n) time. This gives the 0{^/nlog'^ n) query time of their 
data structure. 

We use dense distance graphs in two of our data structures (Sect. [3]and[5|). In both cases it 
is on a variant of recursive decomposition, as discussed above. 

2.3 The Monge Property 

Apxq matrix M satisfies the Monge property if for every two rows i < k and two columns j < i, 
M satisfies Mij -\- M^^ < Mm + M^j . We can find the minimum element of M by transposing, 
negating and reversing M, and using the SMAWK algorithm [1] for row-maxima on the resulting 
totally monotone matrix. If we do not store M explicitly, but are able to retrieve each entry in 
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Figure 1: The distances from X to y Figure 2: A falling staircase matrix 
satisfy the Monge property (left) and an inverse falling staircase 

matrix (right) 



0(1) time this takes 0{p + q) time. Note that if we add a constant to an entire row or to an 
entire column of a matrix with the Monge property, then the property remains. 

Consider two disjoint sets X and Y of consecutive boundary vertices on a boundary walk 
of some piece B. Rank the vertices of X from xi to x^x\ according to their order around the 
boundary walk, and rank the vertices of Y from yi to y^Y\ according to their order in the opposite 
direction around the boundary walk. For i < k and j < £, the shortest path from Xi to yg inside 
B and the shortest path from to yj inside B must cross each other. Let whe a vertex common 
to both paths. Then, dsixi^yj) + dB{xk,ye) < dsixi^w) +dB{w,yj) + dB{xk,w) + dB{w,yt) < 
dsixi, yi) + dB{xk, yj) (see Fig. [1]). Therefore, the matrix M such that Mij = dB{xi,yj) has the 
Monge property. The Monge property was first used explicitly for distance queries in planar 
graphs by [H]. 

A partial matrix is a matrix that may have some blank entries. In a falling staircase matrix 
the non-blank entries are consecutive in each row starting not before the first non-blank entry of 
the previous row and ending at the end of the row (see Fig. [2]), inverse falling staircase matrix 
is defined similarly by exchanging the positions of the non-blanks and the blanks. Aggarwal 
and Klawe [2] find the minimum of an (inverse) falling staircase matrix whose non-blank entries 
satisfy the Monge property in 0{q + p) time by filling the blanks with large enough values to 
create a Monge matrix. 

In Sect.[5]we use this tool for finding the minimum of two staircase matrices whose non-blank 
entries satisfy the Monge property. 

3 Linear- Space Data Structure 

In this section we present a data structure with linear space, almost linear preprocessing time, 
and query time faster than any previous data structure of linear space. We generalize the data 
structure of Fakcharoenphol and Rao |14j by combining recursive decomposition of the graph 
with r-decomposition. This is similar to the way that Mozes and Wulff-Nilsen [28] improved 
the shortest path algorithm of Klein, Mozes and Weimann |23] . Mozes and Sommer |27j have 
independently obtained a similar result. 

We find an r-decomposition of G into p = 0{n/r) pieces, and then we recursively decom- 
pose each piece into p subpieces, until we get to pieces with a single edge. The depth of the 
decomposition is 0(log ?i/ logp) where at level i we have pieces, each of size 0{n/p^) and with 
0{y^njjf) boundary vertices. Constructing this recursive decomposition takes 0(n log n|^|^) 
time. An alternative way to describe this decomposition is to perform a recursive decomposition 
on G while storing only levels klogp for k = 0, . . . , [log n/ log j^J of the recursion tree and the 
leaves of the recursion (the pieces containing single edges) . 

We compute the dense distance graph for the recursive decomposition, in the same way as 
in the data structure of [H] . That is, we compute the distance between every pair of boundary 
vertices in each piece. Using the algorithm of Klein [22] this takes O(relogn) time for each 
level, and a total of 0(n log n|^^) time. The size of dense distance graph over our recursive 
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decomposition is 0{n-^^). 

When a distance query from it to arrives, we use the Dijkstra implementation of |14j to 
answer it. We run the algorithm on the subgraph of the dense distance graph that includes all 
the pieces that contain either u or v, and the p — 1 siblings in the recursive decomposition of 
each such piece. We require the sibling pieces because the shortest path can get out of a piece 
B into a sibling of B without getting out of any piece that contains B. Therefore, the number 
of boundary vertices involved in each distance query is 
Hence the query time using the algorithm of [1^ is 0{py/nlog^ n). 

We conclude that for a planar graph with n vertices and any p G [2,n], we can construct 
in 0(nlog^ n/logp) time a data structure of size 0{n\ogn/ \ogp) that computes the distance 
between any two vertices in 0{p^/n\og^ n) time. If we set p = 2 we get exactly the data 
structure of jl4j . If we set p = for a constant < 5 < e we get: 

Theorem 3.1. For a planar graph with n vertices and any constant e > 0, we can construct 
in 0(n log n) time a data structure of size 0{n) that computes the distance between any two 
vertices in 0{n^^^~^'') time. 

The total time for k distance queries is 0{n log n+kn^^'^'^^) and the required space is 0{n+k). 
This improves the fastest time for k distance queries for k = 0(n^/^~^) and k = a;(logn) 
simultaneously. Among data structures that require only 0{n + k) space, the upper bound on 
the range of k is 0(n^/^~^). 

Fakcharoenphol and Rao [2] noted that Smith suggested that their algorithm can be gen- 
eralized to graphs of bounded genus. If a graph G with bounded vertex degree is embedded in 
an orientable surface of genus g, then [12l|20] showed how to find a planarizing set of 0{^yng) 
edges whose removal from the graph makes the graph planar, in 0{n + g) time. We use the pla- 
narizing set for the first decomposition of the graph, and combine the Dijkstra implementation 
of [H] with standard implementation using a heap for the topmost pieces in the recursion. We 
get that the bounds of Theorem 13.11 apply also to graphs embedded in an orientable surface of 
a fixed genus. 

4 Improved Preprocessing Time for S E [n^/^, n^] 

In this section we present a data structure that matches the space-query time tradeoff of the 
data structures of Djidjev [9l (§5)] and Chen and Xu [8] with the preprocessing time of the data 
structure of Cabello [5]. Our data structure combines parts of the data structures of [9l (§5)] 
and of [S]- 

First, we construct an r-decomposition of G 'mO{n log n) time, for some parameter r 6 (0, n) . 
For each piece B our data structure has three parts: 

(i) The distances dciu^v) and dG{v,u) for every u G dB and v G V{B). 

(ii) A data structure that reports dB{u,v) in 0{y/r) time for u,v £ V{B). 

(iii) For each hole H of B we store dniu, v) for every u G dB[H] and v G V{H) such that v is 
a boundary vertex of some piece contained in H. 

Part (i) is from the data structure of Cabello [5]. The construction of this part requires 
0{n\ogn + r^^'^) time and 0{n + r^/'^) space per piece [3]. Part (ii) was used both by Djidjev [21 
(§5)] and by Cabello [5]. This is the data structure of [3l[9l (§3)] with S = r'^/^, its construction 
takes 0(r^/^) time and space per piece. Part (iii) is from the data structure of [HI (§5)], but 
we construct it more efficiently. We find the distances for this part using the multiple-source 
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Figure 3: The setting in Sect. [5j The vertex u is in the piece and v is in B' . The subgraph J 
is between B and B' . Every path from u to f must contain a vertex from the boundary walks 
X oi B and a vertex from the boundary walk Y of B' . 

shortest paths algorithm of Klein [52] for every boundary walk. The required space per piece 
for part (iii) is 0{n) and the preprocessing time is 0(n log n). 

Since there are 0{n/r) pieces, each with a constant number of holes and 0{^/r) boundary 
vertices, constructing the three parts takes 0((n^/r) log n + riy/r) time and 0{n?/r + Uy/r) 
space. 

Let n, -y be a query pair. We use the data structure of this section to find dG{u,v) in 
0(\/r log r) time. If u and v are in the same piece then we find the distance from u to v using 
parts (i) and (ii) of the data structure with the query algorithm of [5] in 0{y/r) time (see details 
in Sect. 15.21 below). If u and v are in different pieces then we find the distance using parts (i) 
and (iii) with the query algorithm of (§5)] in 0(-y/r log r) time (see details in Appendix IX]) . 

We conclude that for a planar graph with n vertices and any r S (0, n), we can construct in 
0{{n? /r) log n + riy/r) time a data structure of size 0{'n? /r + riy/r) that computes the distance 
between any two vertices in 0{\/r\ogr) time. The sum n^/r + n^/r minimizes at n^/^, and for 
r = -n? / S we get: 

Theorem 4.1. For a planar graph with n vertices and S G [n^/^,n^], we can construct in 
0(5 log n) time a data structure of size 0{S) that computes the distance between any two vertices 
in 0{nlog{n/^/S)/^/S) time. 

5 Improved Query Time for S G [n^/^, n^] 

In this section we present a data structure with an improved query time, for the same range 
of space bounds as in the previous section. In return, the preprocessing time is higher. For 
this purpose we use minimum search in Monge matrices. While previous planar distance data 
structures have taken advantage of the Monge property before, this is the first to use fast 
minimum search in a Monge matrix with the SMAWK algorithm [1]. 

Again, we construct an r-decomposition of G. Assume that we want to find the distance 
from a vertex u to a vertex v that are in two different pieces. Let B and B' be the different 
pieces that contain u and v respectively, let H and H' be the holes of B and B' that contain 
V and u respectively, and let X = dB[H] and Y = dB'[H']. Let J = H D H' he the subgraph 
of G contained both in H and in H'. We assume without lost of generality that J contains the 
infinite face. See Fig. [3l 

The shortest path from u to v must contain a vertex x £ X and a vertex y £ Y {it is 
possible that x = y). We assume that there is no internal vertex of B or B' between x and 
y in this path, since otherwise we can replace x with a later vertex. Therefore, dG{u,v) = 
min^.eX,2/ey{c^G(^i, x) + dj{x, y) + ddy, v)}. 

Our goal then is to find y that minimize daiu., x) + dj{x, y) + dciy, v). For a particular 
order of X and of Y, which we specify below, let M be the matrix such that Mij = dj{xi,yj), 
and N be the matrix such that Nij = dciu, Xi) + dj{xi, yj) + dcivj^v). We show how to order 
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Figure 4: Finding L{x). Only tight edges of J are shown, B and B' are shaded. The leftmost 
path from u' to v' is 6oW, y = i{x), L{x) is its subpath between x and y. Directions of the edges 
are not shown. 

the members of X and Y such that M decomposes into two staircase matrices, each with the 
Monge property. Since dciu^Xi) is fixed for a fixed Xj, and dciyj^v) is fixed for a fixed yj, 
then also consists of two staircase matrices with the Monge property. Thus we can use the 
algorithm of Aggarwal and Klawe [2] to find the minimum entry of A^, which is the desired 
distance. 

For every x £ X we define the leftmost shortest path from x to Y, denoted by L{x) as 
follows. We add to the embedding of J a vertex u' inside B and connect it with an edge to x, 
and a vertex v' inside B' and connect it with edges from every vertex of Y (recall the internal 
vertices of B and B' are not in J). We set the length of all new edges to be 0. An edge 
e = {w, z) G E{J) is called tight if the length of e is equal to dj{u' , z) — dj{u' , w), that is if e is 
on some shortest path from u' to z. We remove all non-tight edges from the graph, and perform 
a left-first search from u' until we find v' (i.e. we perform a depth-first search from u' , and visit 
the edges outgoing from a specific vertex according to their left-to-right order, see also [22]). 
Let L{x) be the path we obtain by removing the first and the last edges of the leftmost path we 
found from u' to v' , and let i{x) S y be the last vertex of L{x). Note that L{x) is a shortest 
path from x to i{x). The reason we added u' is to decide between two paths that diverge at 
X itself, the reason we added v' is to decide between two paths such that one is a prefix of the 
other. Note that L{x) may contain more than one vertex of Y. Moreover, even x £ X OY it 
is not necessarily true that i{x) = x. See Fig. HI 

Fix some arbitrary vertex of X to be xi and rank the other vertices of X in a clockwise 
order. Let y^Y\ = ^{^i)^ and rank the vertices of y in a counterclockwise order. Let P = L{xi). 

Lemma 5.1. Let Xi £ X and yj £Y. There is a shortest path Q in J from Xi to yj such that 
either Q does not cross P, or every prefix of Q crosses P from the left side of P to its right side 
at most once more than it crosses P from its right side to its left side. 

Proof. Assume that every shortest path from Xi to yj in J crosses P, and let Q be such a 
path. First, we can assume that at the first time that Q emanates from P it emanates from 
its right side - if Q emanates from the left side of P at a vertex w, we can replace the suffix 
of P that starts at w with the suffix of Q and get a path from xi to yj which is to the left of 
P, contradicting its definition as L{x) (see Fig. [5|^a)). Second, we may assume that if Q meets 
P at a vertex w, and then again at a vertex w' such that w' is after w also in P, then the 
subpath of Q between w and w' is the same subpath as in P, since otherwise we can replace 
the subpath of Q between w and w' with the subpath of P (see Fig. El^b)). Last we may assume 
that in two consecutive times that Q crosses P it does so from different directions, since if the 
same direction is repeated twice, and by the previous observation the second crossing precedes 
the first crossing in their order in P, then Q must cross itself (see Fig. [5)^c)). From these three 
observations the lemma follows. □ 

Let M' be the partial matrix of M where M^j is non-blank if there is a shortest path in 
J from Xi to yj that does not cross P, and let M" be the partial matrix of M where M-'j is 
non-blank if every shortest path from Xi to yj crosses P. The partial matrix M' has the Monge 
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(a) If Q emanates left from P at w (b) If P and Q contain w and w' in the 
then we can replace the suffix of P same order then we can replace the sub- 
that starts at w with a suffix of Q. path of Q between w and w' with the sub- 



path of P. 





(d) If a shortest path from X to y 
does not cross P [solid lines), or al- 
ways cross P [dashed lines), then M 
has the Monge property (cf. Fig. [I}. 



Figure 5: Crossing of Q and P. 



property, we get this by cutting open J along P and using the claim from Sect. 12.31 (see 
in Fig. [5|^d)). A similar argument (by taking two copies of J open at P and "gluing" the right 
side of one of them to the left side of the other) shows that M" also has the Monge property 
(see in Fig. [5]^d)). The non-blank entries of row i in M' are from l{xi) to y\Y\ = ^{xi), 

and the rest of the row is in M" . The partial matrix M' is a falling staircase matrix and M" 
is an inverse falling staircase matrix, since -L(xj+i), the leftmost path from Xj+i to Y , cannot 
cross the path L(xj) from its right side to its left (this is similar to the illustration in Fig. [5][^a)). 
Let A^' and N" be the corresponding staircase matrices partial to N (with same blank entries 
as M' and M" respectively), both of them have the Monge property. 

In order use the insights above, we take parts (i) and (ii) of the data structure of the previous 
section together with the following two new parts, where X, Y and J are as defined above (part 

(iii) is not necessary): 

(iv) For each two pieces B and B' we store dj{x, y) for every x £ X , y £ Y . 

(v) For each two pieces B and B', and every x £ X, we store i{x). 

Notice that dj{x,y) and i{x) depend on the specific pieces B and B'. Since there are 
n/r pieces, each piece has ^/r boundary vertices, and each vertex is on the boundary of a 
constant number of pieces, the total space required for part (iv) is 0(n^/r) and for part (v) 
is 0(n^/r^/^). This does not increase the total space complexity of the data structure which 
remains 0{'n? /r + riy/r). 

5.1 The Preprocessing Algorithm 

We have an r-decomposition of G obtained by a recursive decomposition, where we decomposed 
each piece into two pieces using a separator and stopped the recursion at pieces of size 0(r), 
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from which we took only the pieces that correspond to leaves of the recursion tree. Now we will 
take all the pieces of the entire recursive decomposition which defined the r-decomposition. We 
build a dense distance graph for G based on this recursive decomposition using the algorithm 
of Fakcharoenphol and Rao [13] with the improvement of Klein [22] in O(nlog^n) time and 
0(n log n) space. 

For two fixed pieces B and B' we use a subgraph of the dense distance graph to compute 
the distances from vertices of X to vertices of Y in J. We should choose carefully a set S of 
pieces of the recursive decomposition to use, we must obey three rules - we cannot take any 
piece of the O(logn) pieces that contain either B or B' , we should cover all the paths from X 
to y, and the total number of boundary vertices of pieces in S should not be too large. 

We start with the entire graph G, which is the root of the recursive decomposition, as the 
single piece in S. As long as there is a piece C in S containing either B or B' in it (C is an 
ancestor of B or B' in the recursive decomposition), we replace C with both of its two children 
in the recursion tree. When we get to B or B' , we remove them from S. At the end of this 
process, S contains O(logn) pieces with 0(-y/n) boundary vertices. All the vertices of X and 
Y are boundary vertices of pieces of S (because otherwise such a vertex will be internal vertex 
of some member of S, which is an ancestor of B or B'), and the paths between them in J 
are covered by S (since S only misses internal vertices of B and B'), as required. Denote the 
subgraph of the dense distance graph that include exactly the pieces of S by D. 

For X G X, we compute part (iv), by computing dj{x,y) for every y GY using the Dijkstra 
implementation of [13] on D in 0{^/nlog^ n) time. 

We use the distances of vertices of D from x to find i{x) for part (v) as well. We add to D a 
vertex u' inside B and connect it to x, and a vertex v' inside B' and connect every vertex of Y 
to it, as described before. Since we added an edge with length from u' to x, for every vertex 
z in D, dj{u',z) = dj{x,z), so we can determine for each edge of D whether it is tight in a 
constant time. Even though D is not planar, we can define a cyclic order on the edges incident 
to each vertex using the embedding of G. The order of edges incident to a specific vertex of D 
is defined by the order of the shortest paths in G that they represent. Due to space constraint 
we give the complete details of this process which takes total time of O(nlog'^n) in Appendix 
[B1 Then, we ignore all non-tight edges in D and find £{x) by finding the leftmost path from u' 
to v' . This takes time proportional to the number of vertices in D, which is 0{^/n). 

We perform the computation of parts (iv) and (v) for every B and B' , and x £ X. For 
every x this computation takes 0{^/nlog^ n) time, and we repeat it O^n^ /r^^"^) times. The 
total time required for constructing the data structure is 0((n^/^/r'^/^) log^ n + riy/r) (note that 
the term nlog^n is dominated by this bound). We note that it is also possible to perform the 
computation in 0((n'^/r^) logn) time using the algorithm of Klein (22] for every pair of pieces, 
however this does not improve the time bound for our range of r. 

5.2 The Query Algorithm 

Let u,v he a query pair. We use the data structure of this section to find dG{u,v) in 0{^/r) 
time. Again, let B, B' he the pieces that contain u, v respectively. 

li B = B' then we can answer the query in 0{y/r) time in same way as in the previous 
data structure, which uses the query algorithm of [5]. Either the shortest path from u to v 
is inside B or the shortest path contains some vertex h G dB. In other words, dG{u,v) = 
inm{dB{u,v),'nimh^QB{dG{u,b) + dcib^v)}}. We retrieve the distance dsiu^v) from part (ii) of 
the data structure in 0{^/r) time. For each h € dB we retrieve dciu, b) + dcib, v) from part (i) 
in 0(1) time. Since \dB\ = 0{y/r), it takes 0{y/r) time to go over all vertices of dB and find 
b E dB that minimizes dG{u,b) + dG{b,v). Then, dG{u,v) is the minimum between dB{u,v) 
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and dciu, b) + dcib, v). 

For the case where u and v are in different pieces, let X and Y be as before. Fix xi to be 
an arbitrary vertex of X and let y^Y\ — K^i)- Let the matrices N, N', N" be as before. We 
compute an entry Nij = dG{u,Xi) + dj{xi,yj) + dciyj^v) in 0(1) time using parts (i) and (iv) 
of the data structure. We determine whether an entry of N is in N' or in N" in 0(1) time using 
part (v). Therefore, we can use the SMAWK algorithm [1] as in [2] to find the minimum value 
in N in 0{\/r) time. This value is the requested distance. 

We conclude that for a planar graph with n vertices and any r G (0, n), we can construct 
in 0((n^/^/r^/^) log^ n + n^/r) time a data structure of size 0(r? jr + n-y/r) that computes the 
distance between any two vertices in 0{^/r) time. As in Sect. [U we set r = ri^ jS and get: 

Theorem 5.2. For a planar graph with n vertices and S G [n^/^,n^], we can construct in 
©((S'^/^/i/n) log^ n) time a data structure of size 0{S) that computes the distance between any 
two vertices in 0{n/\fS) time. 
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Figure 6: Finding 7/3. For 2:3 between xi and X2, the vertex 2/3 is between yi and 1)2, since 
otherwise we get that every path from 3:3 to v (dashed) crosses either the shortest path from 
xi to V or from X2 to at a vertex w. 

A The Query Algorithm of Section [4] when the Vertices Are in 
Two Different Pieces 

Consider the data structure of Sect. |H Let u,v he a query pair, such that u is in a piece B and 
V is in another piece B' . The query algorithm that we describe here is similar to the one of O 
(§5)], and the complete details are given thereH Let H be the hole of B that contains v and let 
H' the hole of B' that contains u. Denote X = dB[H] and Y = dB'[H']. We assume without 
loss of generality that H Ci H' contains the infinite face. 

A shortest path from utov contains some vertex x £ X and some vertex y £Y (it is possible 
that X = y). We may assume that there is no internal vertex of B between x and y in the shortest 
path (since otherwise we can replace x with another vertex of X). Therefore, dciu^v) = 
m.iTix^x,y&{dG{u-,x) + dnix^y) + dG{y,v)}. Next we show how to find minj,(=y{(i//(x, y) + 
dG{y,v)} for every x € X in 0{y/rlogr) time. 

For a fixed vertex Xi £ X it is easy to find yi £ Y that minimizes d}{{xi,yi) + dciyijv) in 
0{^/r) time (the same member of Y may be yt for different members of X), by going over all 
vertices of Y and using parts (iii) and (i) of the data structure. Let yi and 7/2 be the vertices 
that minimizes duixi^yi) + dciyi^v) for xi and X2. There is a shortest path from xi to v that 
contains yi, and similarly a shortest path from X2 to v that contains y2- Let X3 be a vertex 
between xi and X2 in the clockwise order of X starting at xi. There is a vertex 2/3 G ^ that 
minimizes dnix^^ys) + (1/3,7;) located between yi and y2 in the counterclockwise order of 
vertices of Y starting at yi. Since otherwise, every shortest path from X3 to v must cross the 
shortest path from xi or from X2 to v that contains yi or y2, respectively. Assume without loss 
of generality that the shortest path from X3 to v crosses the shortest path from xi to v, and 
let w be the vertex in which the two shortest paths meet. Then, if we replace the suffix of the 
shortest path from xi that begins at w with the suffix of the shortest path from X3 we get a 
shorter path, this is a contradiction (see Fig. [6|). This gives the following algorithm for finding 
yi G Y for every Xi £ X. 

Let xi,X2 be two arbitrary vertices of X. Find yi and 2/2 for xi and X2 by going over all 
vertices of Y. Let x^ be the middle vertex between xi and X2 in the clockwise order of X 
starting at xi. Find 7/3 by going over all vertices of Y between yi and 1/2 in the counterclockwise 
order of vertices of Y starting at yi. Continue recursively for the vertices of X between xi and 
X3 and the vertices of Y between yi and y^, and also for the vertices of X between X3 and X2 
and the vertices of Y between 7/3 and 7/2. Similarly, find yi for every Xi between xi and X2 in 
the counterclockwise order of X starting at xi. 

We conclude that we can find minygy{d//(x, 7/) + dG{y,v)} for every x £ X in 0(-y/r log r) 
time. Now, we go over all vertices of X, and using part (i) of the data structure we find 
minxf^x^yt^Yidciu, x) + dH{x, y) + dG{y-, v)} in 0{^yr) time. The total query time is 0(\/r log r). 

''There is a small difference in our description of the algorithm; instead of Step 1 of [O] (§5)] we use part (i) 
of the data structure which [9] (§5)] does not have. 
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(a) t is the nearest common an- 
cestor of z and z' . 



(b) z is an ancestor of z and 
lies on the boundary of H. 



Figure 7: Two cases where z is to the left of z' . The dotted hnes represent paths in T. 

B A Cyclic Order for Edges Incident to a Vertex of D 

In this appendix we define a cychc order on the edges incident to a specific vertex in the graph 
D, which is a subgraph of the dense distance graph. We use this order in the preprocessing 
algorithm of Sect. 15.11 in order to find i{x) for a boundary vertex x. We define the order of the 
edges such that the leftmost shortest paths from x to y in G, and in D, both end at the same 
vertex of Y {x^Y^lix) and D are defined in Sect. [5]). A vertex u; of D is a boundary vertex of 
more than one piece, however the order between two edges in two different pieces is clear from 
the embedding of G (the pieces of D are pairwise edge disjoint). Therefore, here we define the 
left-to-right order of the edges inside each piece. The left-to-right order of the edges, is in fact 
a left-to-right order of the boundary vertices, because the edges of a piece in the dense distance 
graph connect a vertex on the boundary of the piece to all other vertices on the boundary. We 
define the left-to-right order from w to the other boundary vertices of the piece according to the 
left-to-right order of the leftmost shortest paths from w to the other vertices. This order allows 
us to find l{x) as required. The order that we define does not depend on the specific graph D, 
so we perform the procedure described here only once for every boundary vertex of every piece. 

Let w be a boundary vertex of a piece C . When we compute the distances from w to the 
other vertices of dC for the dense distance graph, we use the algorithm of Klein [22], which 
maintains a dynamic tree [29] that contains the rightmost shortest path from w to every vertex 
of C. Since we are interested in leftmost shortest paths we use a symmetric version of [22], by 
replacing the roles of left and right. Denote this leftmost shortest path tree rooted at w by T . 

Let z and z' be two vertices of dC different from w. We show how to decide in 0(log |C|) 
time which of the two vertices is to the left of the other, with respect to w. Let t be the nearest 
common ancestor of z and z' in T. We can find t and the two edges that lead from it to z and 
to z' in 0(log |C|) time from the dynamic tree [29]. First assume that t ^ z,z' . Consider the 
following three edges incident to t in T - the edge that connects t to its parent i\lt = w then 
we add a dummy edge inside the hole that w lies on its boundary for this purpose), the edge 
that leads from t to z, and the edge that leads from t to z' . The order of these edges around 
t determine the order between z and z' (see Fig. [7)[a)). Now assume without loss of generality 
that t = z. The vertex z lies on the boundary of some hole of C, denote this hole by H. There 
are two edges incident to z on the boundary of H. We can find the two edges when we find the 
piece C. Consider the edge that connects z to its parent in T, the edge that leads from z to z', 
and the place of H among the edges incident to z. If in the clockwise order of edges around z 
starting at the edge that connects z to it parent, the edge that leads to z' is before H , then z' 
is to the left of z, otherwise z is to the left of z' (see Fig. [Zl^b)). 

Since we compare two vertices in 0(log|C|) time, we can use comparison sort to sort the 
0{^J\C\) vertices of dC around the vertex w from left to right in 0{^J\C\\o^ \C\) time. We 
repeat the process for each vertex of dC in a total of 0(|C| log^ |C|) time. For the pieces of a 
single layer of the recursive decomposition, the total time is 0(n log^ n). And for all the pieces 
of the dense distance graph the process takes 0(n log^ n) time. 
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